Pseudo magnetic field in strained graphene: revisited 

M. Ramezani Masir, D. Moldovan, F. M. Peeters 

Departement Fysica, Universiteit Antwerpen 
r/~\ Groenenborgerlaan 1 71, B-2020 Antwerpen, Belgium. 

o 

;-H 

Q^ 

<^ Abstract 

\C) We revisit the theory of the pseudo magnetic field as induced by strain in graphene using the tight- 

^N binding approach. A systematic expansion of the hopping parameter and the deformation of the 

i__i lattice vectors is presented from which we obtain an expression for the pseudo magnetic field for 

'~~| low energy electrons. We generalize and discuss previous results and propose a novel effective 

Cd Hamiltonian. The contributions of the different terms to the pseudo magnetic field expression is 

investigated for a model triaxial strain profile and are compared with the full solution. Our work 
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C/2 suggests that the previous proposed pseudo magnetic field expression is valid up to reasonably 
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high strain (15%) and there is no K-dependent pseudo-magnetic field. 
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I 1. Introduction 

5 Graphene has triggered a broad activity both in fundamental and applied physics and chem- 

istry. The most intriguing feature of this system is their similarity to ultrarelativistic electrons 
and positrons obeying the Dirac equation |[Tl|2l. An interesting prediction is that a geometri- 
cal deformation of the graphene lattice results in a local strain that acts as a pseudo-magnetic 
field on the electronic degrees of freedom and which leads to a pseudo-quantum Hall effect O . 
^v Graphene can sustain very high, up to 25%, elastic strains 101 which leads to a shift in the posi- 

psj tion of the Dirac cones f5l. Deformation due to elastic strain changes the hopping amplitude of 

\^ the carbon atoms and induces an effective vector potential that shifts the Dirac point |7|. With a 

^^ proper geometrical deformation it is possible to create large pseudo-magnetic fields of different 

■^ shapes |l3][8l|9l. It has been predicted that applying strain with triangular symmetry results in an 

^^ uniform pseudo-magnetic field of the order of lOT [6J. Recently it was reported experimentally 

^^ IfTTI that nanobubbles grown on a Pt(lll) surface induce pseudo-magnetic fields of more than 

, . 300 T. Landau quantization of the electronic spectrum was observed by scanning tunneling mi- 

^ croscopy. Thus strain engineering has become a new way to control the electronic properties of 

'"' graphene in [TOl. 

The effective vector potential induced by strain was derived in Refs. flS^T?] and was based on 

C^ a tight binding approach with the important approximation that the local strain does not alter the 

lattice vectors. Very recently it was shown that including the deformation of the lattice vectors 

leads to an extra term for the effective magnetic field which is of the same order of magnitude 

and which differ in the different K points [151 . But later it was shown that this extra term in the 

Preprint submitted to Solid state communication April 29, 2013 



effective vector potential does not have any contribution to the induced pseudo magnetic field 
I.16J and that subsequently there is no different in the different K-points. Furthermore, in Ref. 
ifTTlfTSl it is shown that in the presence of strain the Fermi velocity becomes spatial dependent. 
In this manuscript we revisit the problem and present a systematic study of the different cor- 
rections to the vector potential and compare them with the numerically obtained full pseudo mag- 
netic field. We present the effective Hamiltonian that includes different contributions of strain. 
The previous result for the vector potential and the Fermi velocity reobtained in our systematic 
expansion. As an example we present explicit analytical results for strained graphene as induced 
by a uniaxial and triaxial strain. We find the magnetic field induced by the in-plane deformation 
and compared the different terms for the vector potential 1131 fT4l [TSl with the exact numerical 
results for the pseudo magnetic field. 



2. Strain Field 



The Hamiltonian in the tight-binding approximation considering only the first nearest neighbor 
is given by: 

H - -Y'toalbj + h.c, (1) 



ij 



where to is the hopping parameter and a,- and a- (b, and b.) are the annihilation and creation 
operator for an electron on sublattice A (B). In the presence of lattice deformation the hopping 
parameter t changes due to the changing interatomic distance. The modification of the hopping 
parameter due to strain is given by |5 1, 



t„ = toe 



-m./a-i) 



(2) 



where a is the unstrained nearest neighbor distance, /3 ~ 2 - 3.37 and d„ is the length of the 
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Figure 1: (a) A monolayer graphene lattice, ai = (3/2, V3/2)a and 32 = (3/2, - V3/2)a are the basis vectors and the 
sublattices are connected by (^i = (1/2, yi/2)a,62 = (1/2, - V3/2)fl and(53 = (-l,0)a. (b) The normal (solid black line) 
and deformed (yellow dashed line for full solution and blue dotted for approximated pseudo magnetic field) Brillouin 



zone with K points given by Kj = -^ (0, 1), K2 = -^ f V3, - 1) and K 

3 V3o 3 v3(7 ^ / 
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Strained lattice vector. Using the Fourier transform of the creation and annihilation operators we 
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obtain the strained Hamiltonian as: 



(3) 



dn = (/ + u)6„, where / is the unity matrix and u is the strain tensor. The strain elements of the 
tensor are given by IIT8lfT6l . u - e + co, and it consists of two parts: the Hnear part of strain tensor 
given by, 

1 ( dui duj 



'■' 2 Idx: dxi 



and the rotational part w given by 



1 ( dui duj 



'■' 2 1 dxj dxj 



(4) 



(5) 



On the other hand we can obtain the change in the lattice vectors size as 

dt' — dXjdx'j — dxjdx' + Iduidxj + duidui 

duj dui dui 

I dxidxii H ( 

dxk dxk dxk 



2 r. dui dui dui 
■' di + 2- — dxjdx]^ H dxidxk 



= df +2- 



1 I dui duk\ 1 I dui duk . 
--— + -— dxidxk + :r 3 -J— dxidxk 

2 \ dXk dXi I 2\ dXk dxi 



dui dui 

+ ; — dXidXk 

dxk dxk 



(6) 



2 /_ _ \ du[ du[ 

= d{ +2 £ijdxidxk + coijdxidxk H dxjdxk- 

^ ■'I dxk dxk 

The first symmetric term denoted by 1 correspond with the linear part of the strain and the second 
term denoted by 2 correspond with rotational tensor and has zero contribution to the nearest 
neighbor vector sizes. The hopping changes with carbon-carbon distance but the rotational tensor 
term does not contribute to it. 

First we drive the effective Hamiltonian by expanding Eq. (|3]l up to the first order in strain 
[151 (The second order terms are included in the subsequent discussion but are not listed in the 
expansion of the Hamiltonian because of those expressions are rather involved), 



H = -J^t„ 



11=] 

3 



«=1 

3 





„i(K+q)d„ 





„-i(K+q)d„ ^ 





Zj'" ^'Ka„ 



-/Ka„ 




(1 + /cr^K ■ ua„)(l + io-yq ■ d„) 



(7) 



^ fo ( 1 - ^a„ ■ u ■ a„U -(o- ■ a„)o-,j (1 + icr-K ■ ua„) 
X (1 + icr^q ■ a„ + ia-,q ■ ua„). 
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The different terms of the effective Hamiltonian are shown in Table. 1. The first term is the 
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famous Dirac-Weyl equation. 



//o = - ^ f{) ( -(o- ■ a„)crJ (io-,q ■ a„) = v'ocr ■ 



n=l 

//o = —ihvocr ■ V 
There are three terms in the first order of strain, 
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(8) 



Hi -2^toi -(o- ■ a„)cr, I (cr,q ■ ua„) = vqo- ■ m ■ p. 



2 
the second term is /^-dependent and is given by 



«=i ^" ' (9) 

1 " 

1 = -invQ\ - - " ■ 



Hi - —ihvo I cr ■ M ■ V H — <T V^ ■ u 



Hi^Y^to i—a„ ■ M ■ a„ j I -(o- ■ a„)o-,| (cr-q ■ a„) 



— o- ■ (2m + Tr(u)I) ■ p, (^0) 



H2 = -iK—^o- ■ 2m ■ V + rr(M)/ ■ V + V ■ m + -V ■ rr(M)/ 

and is the same as the term introduced in Refs. ifTTlfTSl . The third and last term is y6-independent 
and is given by 



i/3 = - ^ fo ( -(o" ■ a„)cr,j (/o-jK ■ Ma„) (/o-,q ■ a„ 



., _ .- ) 

= /-^o--(K-M-w)-p *^^^^ 

v%a { 1 

//3 = h-^a- ■ UK ■ M ■ 6>) ■ V + - V ■ (K ■ M ■ w) 

Here v^ = 3ta/2h, cr = (cTi., cr,,) are the Pauli matrices, and co - {-cr,, cr^). In summary we 
can write the full effective Hamiltonian up to the first order in strain, considering both the /3- 
dependent and jS-independent terms as 

Heff^Ho + Hi+H2 + Hi. (12) 

Now the y6-dependent Fermi velocity is replaced by a tensor. 



Vf 



= y°(/ + ^[2M + rr(M)/]) (13) 



which is space-dependent llTTlfTSl . 

Next we derive the pseudo magnetic field induced by strain. The pseudo-magnetic vector 
potential A^j = A^^ + iAy induced by strain is given by 15J, 



1 3 

v = — yf„e-''^''", (14) 



Table 1: Different terms induced by strain in the expansion of the vector potential. Right column indicates the order of 
these terms in the strain (i. e. 0(u^)) and their effect on the different properties. 



n 


Expansion teims 




1 


g-'K-ft, 





2 


-/q-,J„.-*^" 


Dirac equation 


3 


-/q ■ u ■ <y„e-'K'S" 


vp y6-independent 


4 


-5((5„-M-<J„)e-'KA 


Effective vector potential Aq 


5 


'^_{6„--u-6,diS{-a,de-''^-'' 


Fermi velocity v^ jS-dependent (9(1) 


6 


^((J„-S-(J„)(q-S-<y„)e-'K''" 


Fermi velocity vp /^-dependent (9(2) 


7 


-/(K-M-5„)e-'K'S" 


Effective A /^-independent (9(1) 


8 


f (5„ ■«■*„) (K-M-(J„) e-**" 


Effective A y8-dependent 0(2) 


9 


|i(5„-S-5„)2e-'KA 


Effective A yS-dependent 0(2) 


10 


-i(K-M-<J„)V'K*^ 


Effective A y8-independent (9(2) 


11 


-(K-S-(5„)(q.5„)e-'K-«" 


Vf y6-independent (9(1) 


12 


-(K-S-(J„)(q-S-<y„)e-'K'i" 


Vf y6-independent (9(2) 


13 


1 (($„■«■($„) (K-M-5„)(q ■5„K''^''" 


v/r yS-dependent (9(2) 


14 


-|^(5„-M-5„)2(q-<J„K*''" 


v/r yS-dependent (9(2) 



where Vf is the Fermi velocity, t„ are the strained nearest-neighbor hopping parameters. Note 
that Aps is imaginary because strain breaks inversion symmetry in the nearest neighbor hopping. 
The effective pseudo magnetic field induced by strain will shift the K-points as K„ — » K„ + A„ 
(see Fig. [TJb)). Writing the wave vector k with respect to the Dirac cone using k = K + q and 
expanding the exponent and hopping parameter f„ up to second order we find: 



ro 1 1 -H 6t„ + -5A 1 - /K ■ M ■ 5„ - -(K ■ M ■ d„f 



(15) 



K-(5„ 



X (1 - /q ■ 5„ - /q ■ M ■ 6n)e " 



where 6t - -^_6n-u-5n- The effective vector potential is given by q independent terms. Keeping 
the hopping parameters up to second order and expanding e '^■''" we find 



t„e 



/Kd„ 



f < 1 + 6t„ -iK- u ■ 6„ 
[ 1 2 

1 9 1 7 

i6t„K -u-Sn- -(K ■ u ■ 6„Y + -6tf, 



(16) 



„-'K«„ 



The first correction term is the one obtained in Refs. lfT3lfT4l and the second correction term was 
recently added by Kit et al. fTSl . The third term is the new higher order correction term which 
we will add. 
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Considering only the first term we take constant lattice vectors d„ - S„. Using 
the three nearest neighbors vectors in real space (as shown in Fig. ^ 6, =j(l, V3), 

62 - -(1, - vS), S3 - fl(-l,0) and the position of the K-points are given by Ki — — —- (0, 1), 

2 3 VSa 

K2 = — -=- ( Vs, -1 j and K3 = — —- (- Vs, -1 j we obtain the vector potential in terms of the 

3 y3a 3 y3a 

strain tensor elements. 



A, = 



0OyS I Ujcx - Uy 



Ana 



2m.,v 



(17) 



where (pQ - h/e is the flux quantum. 

In order to obtain the correction given in Ref fiS) we need to include the change of the lattice 
vectors with deformation as <!„ = (/+ u)6„. Including this correction we find the following extra 
term to the vector potential for the different K-points, 



\Ki _ 00 4 
2 - 2a3V3 



My,, 



•at j 



\K2 _ 00 

2 2fl I 2 



2m - ^M 
3 ^xy 9 ">'>' 



^2 ~ 2fl 



3 ^xx 



iUx 



Ux 



2V3 

9 ^xy 

2V3, 



(18) 



3 *^xy 

_2„ _ 2V3 

3 ^xx 9 ^xy } 



It is possible to show that this eff'ective vector potential has the form of "^x- We start with 

3foa 

A2 = -- 



{K ■ M ■ s) 



'.J 



3foa 
3r()fl 



l^'^i'^)-" 



V(K ■ u) ■ s 



(19) 



(20) 



and the two components of the vector potential is given by the real and complex part of A2 as 

A.v oc 5 V (K ■ u) 

Ay oc Cjy (K • U) 

and the magnetic field is given by B2 = V x A2 = 0, which shows that there is no K-dependent 
pseudo-magnetic fields. 

Next we include the second order strain part and try to find effective vector potential, 

1 



-i6t„K-u-6„ + -dtl - |(K ■ u ■ 6„f, 



(21) 



h 



h 



and the corresponding vector potential for the diiferent K-points is given by 






JLt J. ^ 



jK2 _ "Ao /? 



'1 



2fl3V3 



-2m^,, + 3 ^[2lUxyUyy + 3 y/3UxxUxy - 3M^^, - UxxUyy 
-3 y/3ul^ + 3UxxUxy - ^/3UyyUxx - 2 ^/3uly + 3UyyUxy 

2ul^. + 3 y/3UxyUyy + 3 yf3UxxUxy + 3u^.y + UxxUyy 
2a 3 V3 (^ -3 V3Mj^ - 3UxxUxy - ^UyyUxx - 2 V3m^^, - 3UyyUxy ^ 

The correction corresponding to 6t^/2 is given by 



/f = 



^f' 



2a 



'0 §;; 



(5mJj - 2m.„m,3, - 4mJ - 3m2 ) 



3^ 
■ 2;r 



UxviUxx + 3Myy) 



(22) 



(23) 



and the vector potential resulting from the last contribution -j-(K-u-6„) in the different K-points 
is given by, 



'3 ~ 3a 9 



2 2 

'^xy ~ ^y}- 

^UxyUyy 






3fl 9 



<Po An 
3fl 9 



3mJ, - 2uly - u^.y + 2 ^l3uxy{uyy - Uxx) 

-2{Uxy- ^/3Uxx)(Uyy- V3m„) 

3m^j - 2u]y - u^.y + 2 yl3uxy(uxx - Uyy) 

-2(Uxy + ^/3Uxx)(Uyy + VSW^y) 



(24) 



This correction is of second order in the strain and is thus important for large strains and the 
corresponding effective field is position dependent. The most important term is In = ^df which 
is K-independent and it is possible to show that the two other terms I[ and It, have a non-zero 
contribution to the vector potential but have zero contribution in pseudo-magnetic field. 

3. Fermi velocity for uniaxial strain 

The tight-binding Hamiltonian for an infinite sheet of graphene is given by, 

m)] 



// = 



/*(k) j 



where. 



/(k) = ^^«^'''"'"- 



(25) 



(26) 



Here, t„ is the strained hopping parameter which is given bv llTTl . 



t„ - toe 



-fitij,, 



(27) 



where w„ = Inla^c - 1- Here Iq - -2.8 eV is the unstrained hopping parameter, /„ is the strained 
distance to the nearest neighbor atom n, a^c - 0.142 nm is the unstrained carbon-carbon distance 
and/? - 3.37 is the strained hopping energy modulation factor. The strained nearest-neighbor 
vectors are given by dn = ( 1 H- u)6n ■ 
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Figure 2: Top: Contour plots of v,i,/vf near the Dirac point for: (a) vq unstrained grapliene, (b) iVj uniaxial zigzag strain 
and (c) Vac uniaxial armchair strain. Bottom: Fermi velocity along the cuts where (d) ky = and (e) k^ = 0, for all there 
cases of the velocity: vq, v-z and VqC. The solid black line indicates the traditional continuum limit Fermi velocity Vf. 
The strain intensity is 10%. 

We calculate the energy spectrum ^(k) of a graphene sheet from the tight-binding Hamilto- 
nian. The velocity can then be obtained as v = Vk£(k). We calculate the velocity for three cases: 
1) unstrained graphene, 2) graphene strained in the zigzag (zz) direction and 3) graphene strained 
in the armchair (ac) direction. The results are shown in Fig. |2] Note that we only consider the 
part of the spectrum that is close to the Dirac point where the continuum limit may be applied 
(up to 300 meV). The velocity obtained from tight-biding is scaled by the traditional continuum 
limit Fermi velocity Vf - -^ . In the case of unstrained graphene from Fig. 2[a), the deviation 
of Vtb from vp is generally smaller than 3%. However, moving to strained graphene (see Figs. 
|2jb,c)), the velocity deviates from Vf by as much as 25%. 



4. Pseudo-magnetic field for triaxial strain 

The displacement of triaxial strain is given by u(r) - (Ujc, Uy), 

Ux — 2cxy, 



"{x^ - y^). 



(28) 



where c is a constant. The corresponding strain tensor M,j(r) = djUj is, 



u{r) - c 



-y) 



(29) 
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Figure 3: Contour plots of: (a) the displacement profile |u(r)| of the triaxial strain; (b) strain distribution, (c) pseudo- 
magnetic field calculated using the the full value of the hopping parameter. The constant of the triaxial strain is: c = 
0.015 nm-'. 

The pseudo-magnetic vector potential induced by strain in graphene is given by Eq. ( [T6] l, and 
the pseudo-magnetic field is then found as Bps = V x Apj. The vector potential depends on the 



strained hopping parameter Eq. (27 i, which can be expanded as. 



ti/to = 1 + ft*" + ftp' + ftf* . 



/I) 



(30) 



(31) 



Usually, only the first order term St. is taken. Here, we will evaluate the effect of the inclusion 
of the higher order terms. 

A contour plot of the displacement profile of the triaxial strain is shown in Figs. [3fa). The 
pseudo-magnetic field in Fig. [3fc) is calculated using the full hopping parameter from Eq. ( [27| . 
The pseudo-magnetic field is mostly homogeneous in the center. Away from the center, the mag- 
nitude of the field follows the triangular shape of the displacement with high magnitudes of the 
pseudo-magnetic field corresponding to locations of large displacement. In Fig. HI we plot the 
pseudo-magnetic field for different approximations of the hopping parameter Eq. pO| . The fig- 
ures are shown in pairs, with the top ones presenting the magnitude of the field and the bottom 
ones presenting the difference between the approximate and the full pseudo-magnetic field cal- 
culated without approximations (see Fig. pjd)). Taking the first order approximation, see Figs. 
Hfd) and (f), results in an almost completely homogeneous pseudo-magnetic field, which shows 
large differences compared to the full solution. Taking the second order approximation results in 
a less homogeneous field which, however, shows a circular symmetry instead of the triangular 
shape of the full field. Calculating the pseudo-magnetic field using the third order approximation 
finally shows the same triangular shape as the exact pseudo-magnetic field. Adding the fourth 
order term further improves the accuracy, but the correct shape has already been achieved with 
the third order approximation. 
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Figure 4: Top: Contour plots of the pseudo-magnetic field generated from different approximations of the hopping 
parameter up to: (a) first, (b) second, (c) third and (d) fourth order Bottom: (e,f,g,h) Diff'erence plots between the 
respective approximations (a,b,c,d) and the field calculated using the full value of the hopping parameter, as in Fig.[3ld). 
The parameters of the triaxial strain is the same as in Fig. [3] 



In Fig. |5]we compare the pseudo-magnetic field approximations at three points as function of 
appHed strain. We varied the c parameter of the triaxial strain from to 0.025 nm ^ The first 
point (A) is located in the center where the strain remains very low (below 0.5%) even for high 
values of c. Because of the low strain, all approximations are able to accurately estimate the 
field. 

Next, we considered point B, where the strain reaches up to 25%. Because of the higher 
strain, the different approximations start to diverge, although the differences aren't very large. 
The different approximations diverge above 15% strain but the differences remain small even 
above 20%. However, point C shows more significant differences. The approximations diverge 
at already for 6% strain. At high strain, the first order approximation significantly underestimates 
the field (by as much as 350 T). Adding the second order term actually results in an even larger 
underestimation of the field. Finally, adding the third order term corrects the field magnitude so 
that it is in good agreement with the full solution. 

From these results we see that a correct estimation of the field in point C is more difficult 
than in point B even though the maximum strain is actually lower in point C. This is because the 
pseudo-magnetic field depends not only on the intensity of the strain, but also on the direction. 
The strain in point B is mostly uniaxial, as it lies exactly along one of the three strain directions 
(see Fig. l5ta)). On the other hand, point B feels a strong influence from both of the top strain 
directions, so it is strongly non-uniaxial. 
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Figure 5: (a) Contour plot of the triaxial strain with three test points marked as A, B and C. The arrows indicate the 
strain directions. The parameters of the tiiaxial strain are the same as in Fig. [3] (b,c,d) The pseudo-magnetic field as a 
function of the strain at the three test points: (b) A, (c) B and (d) C. The field is calculated for dift'erent approximations 
of the hopping parameter from first to third order {B\ to B3 corresponds to i5?''' to &fi^), as well as the full solution 
(Bfiiii for i5f'-'""*). (e,f,g) The differences of the pseudo-magnetic field approximations compared to the full expression 
(AB, = B, - Bfuii) at the three test points: (e) A, (f) B and (g) C. In all cases the stain constant c is scaled from to 
0.025 nm"', as shown on the top .v-axis. The resulting strain at the test point is shown on the bottom A:-axis. 



5. Conclusions 



We investigated the pseudo magnetic field generated by strain using the tight-binding approx- 
imation. The hopping parameter and the deformation of the lattice vectors are expanded up to 
second order in the strain. The contribution of the different terms are compared with the full 
numerical solution for the pseudo magnetic field induced by a model triaxial strain. 

For our numerical calculation a triaxial force is used to strain graphene and we obtained the 
pseudo magnetic field resulting from the different contributions resulting from different expan- 
sion terms and compared the results with the full solution. Numerical results for uniaxial strain 
clearly show that with applying strain the Fermi velocity is spatial dependent. We included the 
second order term in strain in the calculation of the pseudo magnetic field and showed that the 
first order strain is reasonably valid up to 15% strain and that the pseudo magnetic field is the 
same in all K-points. 
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